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The Line of Action of the Resultant Pressure in Discontinuous 

Fluid Motion. 

By S. Brodetsky, M.A., Ph.D., University of Leeds. 
(Communicated by Sir George Greenhill, F.B.S. Received June 7, 1922.) 

1. The problem of any barrier in a fluid stream is best attacked by the 
method due to Levi-Civita,* of which useful accounts, with extensions, are 
given by Cisottif and Brillouin.J The resultant pressure for any barrier has 
been given in terms of the constants defining the barrier ; but the calculations 
required to find the line of action of this pressure have not been carried out. 
It is our object to supply this deficiency here. 

The motion is two-dimensional. Let the complex variable z(EE %+i>y) 

define position in any plane perpendicular to the generators of the barrier, 

the x axis being parallel to the direction of the stream at infinity. We 

define 

_ d<p _ dyfr __ 8<£> _ ^ 3i/r 

ox ay ■• oy ex 

where u, v are the velocity components, and $, ty are the velocity potential and 
stream function respectively. Let w Ez<f> + ity and define £, O, r, so that 

K EE re i$ = — - ; 12 = log t = log r + l0. (1) 

aw "- 

The figure shows the z, vj, f, O planes for such a problem, with the boundary 
conditions. C is the point of bifurcation of the stream-line IC, CA, CA' are 
the stream-lines in contact with the barrier, AJ, A'J' are the free stream- 
lines. We take the standard dimensions to be unit velocity at infinity 
unit value of K^CA + v/CA') as measured in the w plane, and unit density. 
For other values of these quantities we multiply lengths in the z plane by the 
square of ^(^CA + y/CA') in the w plane, and we multiply the resultant 
pressure per unit length of the barrier by the square of the velocity at 
infinity, by the density, and by the square of Ky/CA + x/CA 7 ) i n tne 
w plane. 

; We introduce the transformation in which the variable r = pe i<r is given by 

a/^ = 2"(t— -j— sin<ro, (2) 

where the point G is given by r = € t(r «. It is immediately verified that the 

* 'Rend, del Circ. di Palermo,' vol. 23 (November, 1906). 

t ' Idromeccanica Piana,' vol. 2, Milano (1922). 

% ' Ann. de Chim. et Pliys., 5 vol. 23, pp. 145-230 (1911). 
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barrier AC A 7 is the semicircle p = 1, •— ?r/2 < or < -f 7r/2 in the r plane, while 
the free stream-lines AJ, A' J 7 become the two halves of the diameter along 
the imaginary axis, namely t ■=> ip> — 1 </o < +1. 
If the barrier is plane we have 



XI = 1 Off ( L £ = - ' 



'to 



e 



— KTfjrpJ' 



g—lCTQrj. 



(3) 



If it consists of two planes meeting at an angle \ir at C (measured away 
from the streaming fluid), then 

/ 1 + 6 t<r »T ^ 



In the general problem we write 

/ l + e to "0T 



^ 1 + g t<r °T \X 

\l-e-vror) 



(4) 



O = X log 



+ A x t + -kA 2 T 2 -f f r A 3 T 3 + i^A 4 T 4 + . . . : 



(5) 



= Cii T -f- 1- Afl 2 T 2 + |« 3 T 3 + 1 10,4?* + . . . , (6 ) 

where the A's, «'s are all real, and the expansion in (5) is convergent for 
\t\ < 1. The angle between the tangents at C in the z plane is again A/7jy 
while the bisector of this angle makes with the x axis the angle (marked cr ) :— 

X&o -f Ai sin a + J A 2 cos 2cr 4- -J A3 sin 3cr 4- -4-A4 cos 4cr + . . . . 

If then we have any set of constants we can find the barrier and its. 
position. 
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2. In our notation the components of pressure on unit length of the barrier, 
ll x , By, and the moment M about the origin in the z plane, are given by 

dz 

¥' 



E y + iR x = I 

B'EB 



1 

2 



(dx—icly\ 

B'EB 



M = real part of \ £ 



zcl 



1 

4 



[ (dx 4- idy) \ (dx — uly) 



B 



B' 



^ 



B'EB 



(7) 



v 



the integration referring to any contour like B'EB, where B' is on A' J', B is 
on A J, and the contour lies wholly within the moving fluid. If we express 
the integrals in terms of r we can take for our contour any curve like B'EB 
(see r plane in the figure). It is well known that 



7T 



7T 



E x = — a^, By = - (a 2 + 4ai sin or ), 



(3) 



(see Cisotti, p. 175 ; Brillouin, p. 195). 

As regards the moment Cisotti does not go beyond giving a formula like (7), 
while Brillouin makes a few general remarks : neither works out the explicit 
value of M, and neither defines the meaning of the term " origin " in this 
connection. 

Since the solution (5, 6) gives us really the value of dzjdr, and we have to 
integrate in order to get z in terms of t, it is clear that no origin has yet been 
defined in the z plane. We shall define the origin as follows. Let dz/dr be 
expanded in terms of r ; we shall soon see that the expansion contains terms 
ia-T" 8 , r~ 2 , r _1 , t°, t 1 , t 2 .... Integrate the series without introducing any 
arbitrary constant, so that we get terms in t~ 2 , r"" 1 , log r, t 1 , t 2 , t 8 .... We 
take the origin to correspond to the value of r that makes the expansion thus 
obtained vanish. 

3. If we use the fact that dz/dw = £ and let the symbol ( )- denote 

that in the expression enclosed by the brackets, — c is to be substituted for 
+ 1, we can write M in the form 



M = real part of 



1 dw j 

B'EB 



ctr 



'£* 



<, diu 7 

4 — - dr 

dr 



B 
B' 



(9) 



But 



£ — exp. (ct\T -f I ta 2 7 2 H- ^agr s -f . . .)> 
so that l/£ is obtained from 5 by changing the signs of all the a's. Also 

cl J£_ = l±I 2 (l- h 2^sinc7 .'r---T 2 ) = l±lf (l+e-o T ) (l-e-'^y 



(10) 
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Hence it follows that ^dw/dr begins with a term t~ 3 . Further, if we use 
the notation 

£ — = — 2 (^-2 + ^-2) _ (g-i + tff-i) + (go + *fto) 

6^7 T 3 T 2 ■ T 

+ (ai + *)Si)+ 2(a 3 +^ 3 )T + ..., (11) 

we at onee see that ^dw/dr is got by changing the signs of /3, a alter- 
nately, i.e., 

1 ^ == _ 2(q- 2 — tft-2) , (ai — tiS-i) , (^0 — t/3p) 

— - (ai — fy8i) + 2(a 2 — t^2)r— - .... (12) 

Again, by the definition of the origin, (11) gives 

y^Ldr = ( a ~3-H/3-~2) + (a-i + tff-i ) 
^ rfr t 2 r 

+ (a + ifio) log r ^-(ai + tyS^T + (a 3 + £/3 2 )t 2 + .... (13) 

P 2 P 

-f (a + ^0) (log p + to) + (ai + tfii) pc l(T + («a + iftz) pV 10 " + . . ., (14) 
Hence we deduce that 

£ *? ^ T \ — (a-2 — ^-2) C 2 ttr + (<*!-- tft-i) eW 

^T '/- p 2 ' p 

+ (ao — tj8o)(logp — wr) + (ai — t$i)pe~ L(r -f (a 2 — 1^2) p 2 e- 2t<T + .... (15) 
If we form the product of (14), (15) we get terms in p~ 4 , p~ 3 , p~ 2 , p* 1 , log p, 

P j P 3 P ; • • *• 

Now let us take the contour B'EB in the r plane to be a semicircle whose 
radius p is very small, so that p -* 0. We can then omit all terms with 
positive powers of p. Also the negative powers of p and the log p term 
must necessarily cut out with the corresponding terms in the other part of M. 

Thus in the product ( f — dr ) I f — dr ) we need only take the terms 

\ dr / \ J dr /- 

independent of p, i.e., 

2(a 2 a~2 + $2/3-2) COS 4<r + 2(a2i8-2 — ^-2)82) sill 4<7 + 2(aia~i -f y8i/3-i) COS 2<r 

+ 2 (ai^-i — a-1/81) sin 2o- + (a 2 + /3o 2 ) ^ 2 ; 

and if we put successively o- = tt/2, o- = —tt/2 and subtract, we get zero. 

To get the other part of M we have to multiply (12) into (13), and then 
integrate. We get terms in t~ 5 , t~ 4 , t~~ 3 , t~ 2 , t" 1 , t°, t 1 , t 2 , ... ; r~ 3 logT, 
T -2 i g T ^ T _i i g T? i g Tj T i g Tj T 2 i g T> ( <s When we integrate we get terms 



Resultant Pressure in Discontinuous Fluid Motion. 365 

111 T~\ T~ 3 , T~ 2 , T- 1 , logT, T 1 , T 2 , ..,; T™ 2 (log T + \), T" 1 (log T -f 1), (log t) 2 , 

T(log t— 1), .... If now we make t = pe t(r and let p be constant and small, so 
that p-»Q, the positive powers of p can be omitted, and the negative powers 
of p as well as the log p parts of the terms log t, (log t) 2 must be cancelled by 
corresponding terms in the part of M already considered. This can in fact be 
verified straightforwardly, but it is not worth while reproducing the rather 
complicated algebra here. The terms we have to consider are thus seen to be 

[4i(a 2 /3~ 2 — a» 2 /32)-*2t(ai/3_i — a-ift)] log r + Oo 2 + /3o 2 ) (log t) 2 , 
and we find that 

M = [{(ot 1 (3~ l --a- 1 /3 1 )--2(a 2 /3-2--«-2f32)} <?-{«? + &*) o*f\ 

— ir/2 
= T [(«l/?-l — «-i£l) — 2 (* 2 /3~ 2 — *-2#2)]. 

The form of M shows that we need find only the first five terms in (11), so 
that we have to expand f only as far as t 4 . If we do this and substitute in 
M we set 



to' 
7T 



M = oo { 3«i%2 — «4 -f 4 sin 0-0 (<% 3 -f 2% — 2ccs) -f 1 6& 2 sin V } . . . (16) 

where, in terms of the A J s we have 

cci = Ai + 2\ cos Co, cc 2 = A 2 ~ 2X sin 2o- , 

& 3 = A 3 + 2X cos 3cr , a 4 = A 4 — 2\ sin 4o- , .... 

4. If the barrier is plane then 

A = 1, A± == A 2 = ... = ' } 

so that by (8), (16) we have 

IT 

E x = 7r cos 2 <t , E^ = 7r sin <j cos <to, M = — -j- sin cr cos <r . 

The resultant pressure is therefore 7r cos cr perpendicular to the plane, audita 

distance from the origin is 

I sin or 

in such a sense that the moment is negative, i.e., above the origin in our 
figure when <r is positive. We can find the origin quite easily. By (3), (10) 

dz 1 -f T 2 

and the origin is given by that value of t for which 

\\e w «T* — -~2~ ) + ( T ) + cos o-o . log r = 0, 

so that p = 1 and a is given by the equation 

J sin (2o- + cr )+ 2 sin cr-fcr cos cr == 0. 

The origin is thus on the barrier when it is plane. 
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It is not convenient to solve the equation for a, and in fact this is not 
necessary. The general value of z in terms of t is 

ze~ t<r » = -U e^or 2 — -'-- ■■;■ ■• -f- t — - ) 4- cos <r • log t. 
Hence for the points A, A 7 ; at which r = 4- t, t = — - 1 respectively, we have 
: >! A — ie l(r o 1 2 — i sin <x + ~ cos cr 1 ? £ A / = — ic L<r » 1 2 + J sin cr -}- ^ cos cr J. 

It follows that the width of the barrier is 

4 + 7T COS CTQ, 

while the distance of the origin from the mid-line of the barrier is 

-J- sin o- 
towards A. Hence the resultant pressure meets the barrier at distance 

I sin ctq 

from the mid-line measured towards A. If then the width of the plane 
barrier is I, the distance of the centre of pressure from the mid-line is 

3/ S ^ n a() 

4 " 4 4- IT COS (To 

the well-known result. 

It is of interest to note that the ratio of the distances of the origin and of 
the centre of pressure from the mid- line is 2 : 3 for all values of a . When 
ao = 0, the plane barrier being normal to the stream, the two points coincide 
on the mid-line. 

5. In more general cases the expansion of ^chv/dr in terms of r is often 
impracticable, and we proceed as follows. Using the notation of (11) we see 
that 

f ydtv 2(a-2 + i/3-2) (a-i4-^-i) (a + tfi )\ 7 
1 b-r t ^ ^ o 4- — o — — - far 



vanishes when r = 0. Hence its value for any point r can be written 
rffdw , 2(a- 2 4-*y8-. 3 ) (*-i + iP-i) («o 4- i/3o) ' 

Thus the value of z in terms of r without any arbitrary constants is 
, = ^±ltA 4- (gJzi + ^-i) + (ao + ^ ) logr 



PJ V— ; 4- 2 ( a ^2 4-^- 3 ) , ( a-i-Hft-i) _ («o 4- *ffo) l ^ T /jitn 



Now the integrand in the integral in (17) has no singularities inside the 
circle \r\ =1, and on the circumference of this circle the only possible 
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singularity is at C, i.e., r = e icr * , when X > 1. Hence if t is any point inside 
or on the circumference of this circle, except 0, we can choose for the path of 
integration in (17) any contour lying within the semicircle defined in the 
r plane, and joining the origin r = to the point t. Let us choose the 
straight line lying between the origin r = and the point t. If r is a point 
on the circumference of jr| = 1, bo that p = 1, r = £ t<r , we get 

z = (a- 2 + t/3-. 2 )e- 2t °--f (a-x + ^-i) e~^ + /'(«o -f ^ ) <r 

Jo LV drj 7 ^ pe ia- p° p 2 p J 

where # is now a point on the barrier ACA' in the z plane. 

Let a = ±tt/2, corresponding to the points A, A 7 respectively, and use the 
value of dw/dr given by (10) : we find 

- 1 r 1-D 2 



+ l 1 ^ :: VF(l-~2psino-o + p 2 )?(^) 
Jo L J P 

2(a- 2 +^-2) __ t (a-i + tjS-i ) («o + ^o)' 
P 3 P* p 



*?/>, 



(a- 2 + ^-2) + A («-l + ifi-lj — T) T? («0 + t/3o) 



Jo - J r 

^ — + -s __ jdp, 

where f( v ), ?(-^) are the values of f for T = ±*/> respectively. The 
straight line joining AA' in the z plane can be appropriately called the chord. 
If B is the mid-point of AA', defining the mid-line of the barrier, we have 
at once 

3D = —(«-2-H/8_2) 

ri r 1 — p 2 
+ Jot" " V K1 ~~ 2p sin cru+ ^) f M+C 1 + 2 p sin ^o+p 2 ) r(-^p)] 

_ 2(a- 2 +*ff- 2 )_ fa + *ffo )l rf 
P 3 f> J P " 

A little algebra shows that 

*to = i - I ^ "47?" [(1 — 2/> sin <7 + p 2 ) £ (ip) + (l + 2p sin o- + p 2 ) £( - v )] 

2/o 8+ 4^~ J ^•■* (18) 

and that the value of the integrand in the integral of (18) is zero when /> = 0. 
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If the function £(t) can be conveniently divided into an even function of 
t and an odd function of t, let these be Z (t), Zi (r) respectively. Then 



Z ( — £p) — Z (t/o), 
and equation (18) takes the form 



Z x ( — ip) =— Zi(tp), 



^D = 



4 







— : — JL. slno-o . Zi (tp)- 



1 , »i a + * O2 + 4% sin cr )\ , n qv 
_+ _ jdp. (19) 

The co-ordinates %>, yj> of the point D are given by (18) or (19), whichever 
is more convenient. Thus for the plane barrier we use (3), and we get 

3 D = J (1 — <? 2to ») = -| sin (Tq (sin 0*0 — 1 cos <r ), 

agreeing with the results already obtained. 

In more complicated problems we can find x D , y D by means of numerical or 
graphical integration, remembering that the integrand in the integral of (18) 
or (19) is zero when p = 0. 

6. Again, by the equations of § 5 we have 



H^a-SaO 



z 






1 — p 



^[(l-2psi.ii(ro + p 2 )?(^)-(l + 2psmcT ( )-+-p 2 )?(-^)] 







4p° 

p 2 J 

Hence we find 

I (SA — 3A') = «A — «D = «D — «A' 

= — (sin (j-o — 5 ^1) — 5 («a + 4% sin <r — £% 3 ) 



2 



8 



f 1 fhzEl [(l - 2p sin cro + p 2 l £ (ip) - (1 + 2p sin <r + p 2 ) f ( - *p)] 

Jo L V 



4- sin 00 



* ,, I 



2 



% Mp ? 



(20 



P' 



or 



(sin o-o ~~ o a - 3 (% + 4«! sm <r - t^i-) 

2 o 



1 n-p* « 



j X.-P 



i-p 2 . 

r en- 



sin 0*0 — ~ % 



, Zi Up) — —£ sin cr . Zo (/.p) -f 
Q L 2p* P 



2 >clp... (21) 



where Z 0? Zi have the same meanings as In § 5, and the value of the integrand 
in (20) or (21), when p = 0, Is 

■^V (6«i 2 sin 0-0 + 12 sin cr + 3r/i«-2 + 6ta 2 sin <r — 2«a 3 — w*i 3 )- (22) 
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We can thus find the positions of A, A 7 relatively to D in the z plane, and 
also the origin. The barrier itself can be drawn by means of its curvature in 
terms of direction, and the problem is completely solved. 

7. In some problems it may be of interest to have the distance of the line 

of action of the resultant pressure from the mid-line of the barrier. This 

is found as follows : Using E^, ~R y% M as already defined in § 2, the moment of 

the pressure about D is 

M— xjtf&y + yvT&x, 

in the anti-clockwise sense in our figure. Hence, by (8), (16), the distance of 

the resultant pressure measured from D towards A is 

1^4 + (ag — fli)sin o- — 2a 2 sin 2 cr — j^ 3 sin cr — § a^a 2 + (a 2 + 4gi sin <r ) a? D — fli 2 #p . 

{a^ + (a 2 + 4^i sin cr ) 2 } 1/2 

(23) 

It may also be useful to know that the resultant pressure meets the line 

through D perpendicular to the direction of the stream at distance 

a 2 + 4#i sin cro 



JL ^ > / ^3 ~~ fll 

8 «1 2 \ Ml 2 



^a x ) sin <x — 2 -— 2 sin 2 <7 — 4- &2 + 



ai' 



xv—yv 



2X sin 4ctq, ..... 



(24) 
from D, measured positive on the side of A. 

8. In the case of the problem defined by (4), so that the barrier consists of 

two planes meeting at C at an angle of X77-, whose bisector makes an angle 

of cr with the direction of the stream at infinity, we have 

&i = 2Xeosa- , a 2 = — 2Xsin2cr , a% = 2Xcos3cr , a 4 = — 

Equations (8), (16) give 

Ea; = 7rX 2 COS 2 CT , E y = 7T A, Sin <7 COS <7 , 

M = ^7rX sin cr COS cr {2(1— X 2 )cOS 2 Gr — 1}. 

Further, substituting for £in (18), (20), we find 

— x - •* - ; 3 (1 — 2p sin (T -{-p 2 ) cos (2X tan x y 



x 



p COS (To 



+ (1 + 2p sin o-o + P 2 ) cos ( 2X tan" 1 ^l 00 ^ 

\ 1 + psm 

X 2 COS 2 (To 



+ p sin cr /J 2p 3 



p sm cr 0/ 
1 



+ 



y» 



p 



dp, 



o 



4/r 



(1 — 2p sin cro -f- p 2 ) sin ( 2X tan*" 1 



p COS cr 



1 



•p sm go/ 



(l + 2psina- Q + p 2 )sinf2Xtan- 1 - P 008 ,^ 



V 



1 + p sin cro/ 



. X sin (To cos o-o! 7 

-r ' >■ Ctp, 

P J 
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\ (^a —%a>) = — sin <to — \ ^ sm a o cos ctq 



p' 







4p 



3 



(1 — 2/} sin cr -j-p 2 )cos ( 2A, tan 1 ~r 



p COS <7o 

— -p sin cr 0/ 



i(y A - 



— (1 -f-2p sin o- -f p 2 )cos (2Xtan~ 1 -JL_i_l£L_ 

• \ 1 + psm €r . 

■ y ) z= A, cos o- 4- i^X 2 cos2 °"o 



+ 



sin <xo 

/>2 



d^p, 



ri 



I A „3 



P' 



Jo *- 



4p ; 



/ o cos 

(1 — 2/} sin cro + p 2 ) sin i 2X tan" 1 ~- — -^ 

\ t — p si 



p COS CTq 

p sin cro. 



+ (l + 2psincJo + P 2 )sin (2\tan 1 ~ ^— . 

' \ 1 -f p sin cro/ 



X cos o- 



/o* J 



Vdp, 



the integrands in the integrals of x-q, t/n being zero when p = 0, while those in 
l-(x A —x A >), ^OjA — yA') have the values 

sin <tq, |- (X — X 3 ) cos 3<r — ^ (X + X 3 ) cos <r 

when p — 0, The knowledge of these values of the integrands is useful when 
numerical or graphical integration has to be resorted to. 

Starting with any point in the z plane as D, and with a given direction of 
the x axis, namely, that of the stream at infinity, xjy, y D fix the origin ; 
h'(%A—%A')> i {V^—Va.') nx the points A, A'. The angles that CA, CA' make 
with the positive direction of the x axis are 

We can therefore draw CA, CA', and the geometry of the z plane is com- 
pletely determined. The line of action of the resultant pressure makes 
an angle 



tan " 



„! / tan cr 



X 



with the x axis, and is at distance 



\ sin <r - 



2 (1 — X 2 ) cos 2 <x •— 1 



(sin 2 a + X 2 cos 2 (7o) 1/2 

from the origin, in such a sense that the moment is anti-clockwise when this 
distance is positive. 

9. A curved barrier of some interest is that defined by 

/ 1 4~ e lcr °T \ ( 1 + l\e l<Xl T \-° os °o/cos o-! sin ai 

where £i, ai, cr 0? en are connected by the relations 

ti = tan ^ ai, sin <r = sin a\ sin <n. 

Since the value of X is unity, the barrier is continuous at C. For different 
values of ai, en, we get first approximations to barriers in the form of 



Resultant Pressure in Discontinuous Fluid Motion. 37 1 






lliptic cylinders of various eccentricities and in different orientations 
relatively to the direction of the stream at infinity. We now have 

2 cos <t , 

cii = 2 cos cr — h cos cri, 

cos cri sm ax 

a 2 = —2 sin 2cr 4 ^-^ — £i 2 sin 2cr b . . . 

cos <xi sm cci 

and equations (8), (16) give 

-13 o / COS C6i \ 2 ^ . / COS ai \ 2 
K x = 7T COS 2 (To 7- — J , K v = 7T Sill O" COS CTq — . 

\1 4- cos «i/ \1 + cos ai/ 

COS cii 

M = — I-7T sin cro cos cr -: =-r„ (1 + 8 cos ai 4- cos 2 ai + 4 sin 2 <7 ) ; 

(1-j- COS Oiif 

while the line of action of the resultant pressure is at distance 

— x ^_± — JliL-- (1 + 8 cos cti + cos 2 ^i -f 4 sin 2 ctq) 
(1 4- cos ai) 2 

from the origin, in such a sense that the moment is anti-clockwise when this 
distance is positive. 

The geometry of the z plane is given by the direction at any point of the 
harrier, and the corresponding radius of curvature K, both expressed in 
terms of a as a parameter : — 

6 = + 5 +ff0 £2i£L_ {tan-i ( fr* n fr + ?> ) 

~ 2 COS CTi Sill ai L \l+^l COS (cT + CTl)/ 

„-, / ^ sin (cr— -cri) \ 
4- tan 2/ ] 



\1 — ^iCOS(o"~0"i)/ J ' 

+ according as a ^ cr , and 



p 2(l + COScr + cro) . — 

it = — . (\ J r sin ai C os cr + cri) 



COS (To COS ai 



(1 — sin u\ cos cr — ci) 



1 + sin ai cos cr 4- <7i\- C0S °-«/ 2cos ^ sina \ 



1 — sin «i cos <j — - cri/ 
We also have 



■#d = i — 



ri r i_ p 2 

4p 3 



[(1 — 2 p sin cr + p 2 ) cos 0i 4- (1 + 2 p sin <x 4- p 2 ) cos 2 ] 



1 cos 2 - / ™° — X2 

"~"7; — ?,4~ 



2p 3 ' p \l4~cos ai/ J 

y D = — 1 - T ^-[(l-~-2psina-o4-p 2 )sin0i--(l4-2psino-o + p 2 )sin02] 
Jo L 4 P 



sin cr Q cos cr / cos cci \ 2 , , 
p \l4-cosai/ ( 
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where 

6 = 2 tan- 1 p G0S ^ — -^21£iL_ tan- 1 f-A£i^£I-J\ 

1 — -p sin (7d cob <7i sin « x '" \ 1 — £ip sin a-\ /' 

2 = 2 tan- 1 _££2LfJL-._ 2 GQS ^ tan -~i. / ^ip cos en % 

1 4- jo si n (7 cos o*i sin ai * \1 -j- ti p sin ai- 

ctXXU. 

2 (^A—^aO = — Sill O- (l +~ COS <r ) _— ±- 

^ \l + cos «!,• 
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j "7~T"f-^ ~" ^i° s * n ^ + p 2 ) cos ^i "~(1 + 2p sin cr + p 2 )cos # 2 ] "i L m 
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i(?/A-2/V) = cos a (1 4- ™ cos <r ) 



2 \l~f-cos^i, 

~~ 1 "T3 [ (^ ~~ ^/> sin o- + p 2 ) sin ft + ( 1 + 2 p sin cj ~f p 2 ) sin O2] 

Jo L ^p 
_co ! o_o/' cos^yi 

p 2 \l+COSai/J 

These give the position of the origin referred to the point I), and also the 
distance AA' as a verification of the barrier itself, when plotted by means of 
the values of R, 6 in terms of o\ 
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